Entangling two defects via the surrounding crystal 
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We theoretically predict that two defects embedded in a crystalline structure, a harmonic chain, 
can be entangled through coupling with the crystal. Entanglement is found at long times for suffi- 
ciently cold chains and for a certain class of initial, separable states of the defects. The calculation 
is performed for a specific physical realization based on trapped ions, but can be extended to other 
systems exhibiting spatial order, such as optical lattices or cavity arrays in circuit QED. The pre- 
dicted entanglement is robust against fluctuations of the chain parameters, such as its size, and 
of the elapsed time, as long as finite size effects can be neglected. Entanglement results from the 
interplay between localized modes in the regime in which the chain acts as a thermal bath for each 
individual defect, and is an example of quantum reservoir engineering. 

PACS numbers: 03.67.Bg, 03.65.Yz, 42.50.Dv, 03.67.Mn 



Entanglement is a quantum mechanical property, 
which has no classical counterpart. One of its peculiar 
features are the correlations of the measurement outcome 
between two entangled objects, even when they are at 
large distances These properties make entanglement 
a precious resource for quantum communication proto- 
cols and quantum metrology applications @- On the 
other hand, its quantum mechanical nature makes it frag- 
ile against coupling with the external world 0] , especially 
when the spatial distance between the two objects is in- 
creased (j46(- Understanding the dynamics of entangle- 
ment in open systems has recently become an active area 
of research focusing on thermalisation in quantum sys- 
tems , and is awaiting experimental verifications. 

Recently the idea has been put forward that the cou- 
pling with a reservoir can allow the robust preparation 
of a physical object in a nonclassical state and can per- 
mit its quantum coherent manipulation [ll-14|. These 
concepts have been applied in order to create entangled 
states of matter [lij], perform quantum simulations with 
trapped ions [l6| , and they are the basis for protocols for 
quantum networks (l7l . [Tq and quantum metrology flif . 
First analyses of the basic requirements for the onset of 
these dynamics, such as the properties of the system-bath 
coupling, have recently been performed |20l - |26| . Condi- 
tions for observing entanglement in optomechanical and 
biological systems are also being discussed 0, Hfjj . 

In Refs. [Ill H3 some of the current authors investi- 
gated a paradigmatic model in which two defects couple 
to a linear chain of oscillators. Generic conditions were 
identified under which entanglement can be established 
between the defects via the common reservoir, which oth- 
erwise would lead to thermalization of a single defect [2^] . 

In this Letter we build on these ideas and show that 
the capability to tune the frequency of the system (the 
defects) in a non-Markovian reservoir allows one to gen- 
erate entanglement between remote systems, even when 
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FIG. 1: (Color online) Two defects of larger mass are embed- 
ded into a linear chain of ions. A laser standing wave illumi- 
nates the defect ions and couples their transverse and axial 
displacements via the mechanical effects of light. For a certain 
class of initial separable states, the defects transverse motion 
becomes entangled by the coupling with the axial phonons of 
the chain. 



the reservoir acts as a thermal bath in the presence of a 
single defect. Memory effects are due to the presence of 
the second defect, so that entanglement is generated by 
means of localized modes of the chain which involve the 
defects and is found for time scales over which the finite- 
size effects of the bulk can be neglected. Our discussion 
here focuses on a specific example, namely, a chain of 
trapped ions which includes two ions of different species 
(see Fig. [I]). Nevertheless, the dynamics we discuss are of 
general nature and can be encountered in systems such 
as dipolar gases in optical lattices 3^, 31|, optomechan- 
ical systems [24|], and cavity arrays in circuit quantum 
electrodynamics J32j. We also mention that a proposal 
for generating entanglement in a chain of 3 ions using 
laser cooling has been discussed in Ref. [33| . In this Let- 
ter we identify a mechanism, that allows one to generate 
entanglement between two ions in a chain of any size, 
and which is solely generated by the coupling with the 
external environment. In a specific limit, the dynamics 
discussed here reproduce the predictions in [33]. 

We consider a chain of ions of charge Q and mass 
to, that is aligned along the z axis and contains two 
defect ions, as sketched in Fig. [TJ The defects have 
same charge Q but mass M > m and are assumed to 
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be trapped sufficiently far away from the chain edges. 
The chain is stabilized by a linear Paul trap, whose 
effective potential is described by the secular potential 

Wrj) = (U\\zf + UxAtf + Vj))/ 2 - Hcrc > 3 = 1, • • • ,N 
labels the ions along the chain, rj = (xj , yj , Zj ) denotes 
their positions, and U\\ and U_\_,j = (Uo/mj — Uu)/2 deter- 
mine the strength of the axial and transverse potentials. 
Here the mass dependence of U±,j is due to the radio- 
frequency potential which creates the transverse confine- 
ment [341 ]. This property, which is specific to ion traps, 
implies that the defects transverse motion is a localized 
oscillation in the chain for sufficiently large mass ratios 
M/m [H. 

The ions are assumed to have been laser cooled to a 
temperature T such that they perform harmonic vibra- 
tions about their respective equilibrium positions rj -* = 

(0,0, zj ''), which arc determined by the balance of the 
trap force and the Coulomb repulsion. In the harmonic 
approximation, axial and transverse excitations are de- 
coupled and the dynamics is governed by the quadratic 



Hamiltonian Hq = Hh + H 



(x) 



with 
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(1) 



(2) 

Here, pj a is the a = x,y, z component 

of the momentum pj of ion j andJCjj = 2Q 2 /\z^ — z^\ 3 
is the coupling due to the Coulomb repulsion (36|. The 
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Hamiltonian term is found from H^' by replacing 
Xj —> yj. We denote the location of the defect ions by 
j = .7i , J2 and assume that rrij 1 — nij 2 — M. The mutual 
distance between the ions is chosen to be small compared 
to the chain length, ji —j\ <C N and we require 1 -C ji < 

h < N. 

Let us first analyze the spectrum of the Hamilto- 
nian Hq. For simplicity, we assume equally-spaced ax- 
ial equilibrium positions with constant interparticle dis- 
tance a = Zj^ — z^°\ Such a situation can be realised 
in optical lattices and exists in the central region of long 
ion chains [37j and in anharmonic axial potentials [38j ]. 
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Fig. 2(a) displays the spectra of the axial and transverse 
modes for the specific case of M ~ 2.87m, which corre- 
sponds to In + ions embedded in a Ca + chain [39|. We 
find two degenerate normal mode frequencies for each 
transverse spectrum, stemming from the localized vibra- 
tions at the positions of the defects, which are separate 
from the respective transverse branch. If the frequency 
gap between localized modes and the transverse branch 
is sufficiently large, these modes coincide approximately 
with the defects transverse vibrations. In this limit the 



defects transverse modes are strongly coupled with each 
other via the Coulomb interaction, and weakly coupled 
to the rest of the chain (40l |. 

Using the isolated frequencies and the long-range 
Coulomb interaction, entanglement between the two 
transverse modes of the defects can be generated by 
preparing each mode in a squeezed state and letting 
the dynamics evolve: the defects will become entangled 
in the same way a beam splitter entangles two optical 
modes which are prepared in a squeezed state [4l|. We 
assume that the chain is initially in a thermal state at 
temperature T while the initial state of the transverse 
modes are identical separable squeezed vacuum states 



along the x direction with variances Axj t ,- 2 



-,2.s 



and Apj i 



e- 2s /2. Here x = ^/MjU\\ 



/2 



is 



the size of the defect ground state in the axial poten- 
tial, po = H/xq is the associated momentum, and s is 
the real- valued squeezing parameter [4l|. The defects' 
state remains Gaussian under the evolution described by 
the quadratic Hamiltonian, Eq.(H|), and is hence fully 
characterized by the first moments and the co- 

variance matrix = | (£i£j + £j£i) — (&) where 
£ = (x J1 /x ,p XJ1 /p ,x j2 /xo,Px,j 2 /Po) [H- Entangle- 
ment between the defect modes can then be quanti- 
fied by means of the logarithmic negativity, Ejy = 
max{0, — ln(2£>_)}, where i>- is the smallest symplectic 
eigenvalue of the partial transpose of the covariance ma- 
trix E 0, HH. Figure [2(b)] displays Ejyasa function of 
time when the defects are at a distance d = la for initial 
squeezing of s = 0, 0.5. One can see that En grows as a 
function of time and increases as s is increased (see grey 
curves). Surprisingly we also find entanglement for s = 0, 
which can be explained by the fact that the defects are 
not eigenmodes of the transverse spectrum. Numerical 
simulations show that this entanglement decreases as the 
mutual distance between the defects is increased and it 
can therefore be understood as a consequence of the form 
of the transverse eigenmodes. 

We will now demonstrate that entanglement between 
the defects can be generated and substantially enhanced 
by coupling the defect modes to the axial phonons of 
the chain. For this purpose we consider an interaction 
which couples the defects' transverse displacements with 
the corresponding axial displacements. The chain dy- 
namics is then given by the Hamiltonian H = Ho + Hi(t), 
where 



= ^ [(; 



Ji 



(3) 



with 7(£) = 7 0(t) being an effective coupling strength 
and 0(f) the Heaviside function. This coupling could be 
realized, for instance, using a standing-wave laser field in 
the x — z plane, resonant with an optical transition of 
the defect ions, and which has nodes at their equilibrium 
positions. The Lamb-Dicke regime is also assumed. At 
t > a displacement of each defect along x excites a wave 
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FIG. 2: (Color online) Results for two In + ions embedded in a chain of 50 Ca + ions with a mutual distance of d = 7a. The mass 
ratio is M/rn = 2.87 [3y]. (a) (Doubly-degenerate) eigenfrequencies ui^ of H± (solid) and Hi\ (dashed line) as a function of the 
quasimomentum k (in units of n/a). The frequencies are in units of Uii = y/Uii/m. The parameters are U±(rn)/U« — 2000 and 
u ii — 2tt x 26.1kHz. The isolated transverse frequencies (see dashed circle) almost coincide with the frequency of the defects 
transverse potential, (b) Logarithmic negativity as a function of time t, in units oiui^ 1 . Black (grey) curves represent the case 

with laser coupling switched on (off), with coupling strength 7 = 18.2mojy . Labels indicate the initial squeezing of the defects 
vibrations, while the rest of the chain is initially prepared at T — 8&j||/fcfl (here: T = 10/iK). (c) Spectral density, in units of 
mojp as a function of the frequency, in units of u)\\. The solid (dashed) line correspond to J+(oj) (J_(oj)). The arrow indicates 
the transverse trap frequency of the defects chosen for the solid curves in (b) [46| ]. 



packet of quasi-resonant axial phononic modes (note that 
displacements along y are decoupled and will be ignored 
from now on). In the following we numerically evaluate 
the defects' dynamics, starting from the formal solution 
of the coupled Hcisenbcrg equations of transverse, axial, 
and defect modes of Hamiltonian H. 

Let us, without loss of generality, consider a sim- 
plified situation, where the defects are symmetrically 
placed with respect to the trap center and where the 
Hamiltonian ([2]) is therefore invariant under the ex- 
change of the j-th and the (N — j + l)-th ion (j € 
{1, ■ ■ • , N/2}). By performing a change of coordinates, 
qf = {qN-j+i ± Qj)/V2 and X± = (x jl ± x h )/\/2, one 
can immediately see that the defects' transverse center- 
of-mass (COM) displacement, X+, couples to the ax- 
ial displacements ', and the same holds for X- and 
qj . The action of the axial vibrations on the dynam- 
ics of the defects' collective variable can be character- 
ized in terms of the corresponding spectral densities 
J^{u) = Efii n{^) 2 /{2muj±)5(uj - uf), where 7* is 
the coupling strength of the coordinate X± to the fc-th 
normal mode of the corresponding axial reservoir with 
eigenfrequency tu^ = uJk [13, EH]. Fig. 2(c) displays 
the spectral densities J±(u) for d = 7a. The appear- 
ance of frequency values ui^ = ojg , at which either the 
spectral density J + or J_ vanishes, is a signature of the 
coupling between the defects and the excitations of the 
ions between them. This can give rise to localized exci- 



tations when the frequencies of the defects can be tuned 
to coincide with one of the zeros, which we denote by 
\ ~ y/U_]_(M)/M 0. Under these conditions, 
one generates a decoherence free subspace for the relative 
motion of the defects, corresponding to an cigenmode in- 
volving the defects and the ions of the chain localized 
between them. Entanglement between the defects can 
thus be generated for an initial state in which both are 
in a squeezed state while the rest of the chain is in a ther- 
mal state. In fact, the dynamics leads to thermalization 
of the defects' COM motion, X+, and destroys all initial 
correlations between the defects' COM and relative mo- 
tion, while the relative coordinate preserves part of the 
initial squeezing. Sufficiently large initial squeezing and 
low initial temperature can therefore lead to two-mode 
squeezing of the defects' transverse motion [28 , 47 1 , and 
therefore to entanglement [48}. We note that, while this 
argument is based on reflection symmetry in the chain 
about the center, such symmetry is irrelevant as long as 
the distance of the defects from the edges of the chain is 
large and their mutual distance is much smaller than the 
chain size. In this limit, our model is valid for any posi- 
tioning of the defects, for time scales over which finite-size 
effects can be neglected. 



The solid curves in Fig. 2(b) show the entanglement 
created in a chain of N = 50 ions, with a coupling 
strength of 7 = 18.2mw^ and two initial values for the 
squeezing. The curves are shown for times over which 
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finite-size effects can be neglected. Entanglement is built 
up after a transient time and reaches values which are 
about an order of magnitude larger than the values found 
in the absence of the coupling laser. Larger entanglement 
is found again by increasing the initial squeezing and also 
decreasing the temperature of the chain (not shown). Let 
us remark that, while cooling a large chain to ultralow 
temperatures is a challenging task, the basic requirement 
for observing the dynamics predicted here is that the ax- 
ial mode of the chain resonant with the defect frequency 
is prepared in the ground state. In a chain of trapped ions 
this could be realized, for example, by cooling the mo- 
tion of the defect ions while the laser coupling transverse 
and axial displacement is switched on (49(. For systems 
of ultracold gases this should not be a critical require- 
ment, considering that the temperatures are of the order 
of few liK. We also notice that noise will not significantly 
affect the predicted dynamics, provided the chain eigen- 
mode couples only weakly to the environment. If, for 
instance, the zero of the spectral density corresponds to 
an eigenmode with short wave wavelength, it will only 
be weakly heated by processes such as patch potentials 
at the electrodes of ion traps [E(j. Entanglement between 
the defects will therefore be observable as long as the ax- 
ial mode of the chain, corresponding to the zero of the 
spectral density, is cold and protected from external noise 
sources. 

We finally demonstrate that the generation of entan- 
glement is not a finite-size effect, and can indeed be found 
in chains of increasing size. Fig. [3] displays the logarith- 
mic negativity as a function of time for various chain sizes 
and a mutual distance of d = 15a between the defects. 
Here we assume that the COM motion is decoupled by 
tuning its frequency to a value at which the correspond- 
ing spectral density vanishes (see inset). The curves are 
displayed for times that are shorter than the revival time 
and have been rescaled by the size-dependent frequency 
W|| = w re f y/log N /N , where oj ro f is a constant and the 
size dependent factor corresponds to a scaling where the 
interparticle distance is kept constant in the center of a 
linear Paul trap. For this thermodynamic limit, the ex- 
citation spectrum of the chain is well defined |36| . One 
observes that the logaritmic negativity oscillates about 
a (quasi) stationary state, whose value is independent of 
the chain size. For increasing chain sizes the time window 
over which this entanglement is found also increases. 

The resulting entanglement can be measured by ex- 
tending the method developed in Ref. [Hl[ to two 
Gaussian modes. The procedure consists of coupling 
the transverse oscillation of each defect with an an- 
cillary qubit, such as an electronic transition of the 
defect ion. This can be done using the Hamilto- 
nian Hj nt = //;/,;/ irrj i '" 1 + a]e 1 ^*), with j G 

{31,32} an d where dj = (xj/xo + ip X j/po)/ V% anni- 
hilates a phonon of the defect transverse excitation 
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FIG. 3: (Color online) Logarithmic negativity as a function 
of time when d = 15a for a chain with N = 300 (dashed), 
500 (dotted) and 700 (solid line) ions (Note that all three 
curves are lying on top of each other). The chain temperatue 
is T — 0, the defects' initial squeezing is s — 0.5. Time is in 
units of u)7 , where = cj rc f\/logiTV jN and ui le t — 2ty x 659.6 
kHz [3^ ]. The other parameters are M/m = 2.87 and 7 = 
0.043 muifef. Inset: corresponding spectral densities, in units 
of mio\ . 



at frequency Clj = y/U±_{M) /M . The frequencies gj 
denote the Rabi coupling with the internal transition 
and <Jj' v ' z are the Pauli operators. The expectation 
value (T) = (a x (g) a x - a y <g> a y + ia x ® a y + ia y ® a x ), 
gives the characteristic function xC^iiAO of a two- 
mode continuous variable system, where ^2 = 
2i J Q dt'gj 1 j 2 (t')e ia ^^ t ' . The reconstruction of the os- 
cillators state and thus of its entanglement properties 
is hence granted by properly designed coupling profiles. 
Note that, as Gaussian states are fully described by their 
first and second order moments, it is sufficient to probe 
their characteristic function close to the phase-space ori- 
gin. 

In conclusion, we have shown that robust entanglement 
can be created between two systems coupled by an ex- 
ternal environment and identified the basic requirements. 
These include the capability to tune the defect frequency 
across a region of the spectrum, which allows entangle- 
ment to be mediated by an eigenmode of the reservoir, 
which must be sufficiently cold and protected from ex- 
ternal noise sources. We have shown that these predic- 
tions can be tested in an experimental setting, focussing 
on the specific case of trapped ions. Nevertheless, due 
to their generality these dynamics can also be observed 
in other physical systems, such as cavity arrays, circuit 
QED and optomechanical systems. Our study sheds new 
light on the dynamics of thermalization and the role of 
the reservoir in establishing quantum correlations, and 
contributes to applications for scalable quantum techno- 
logical implementations. 
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